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CHAPTER 1 


INTRODUCTION 


The handling or flying qualities of a piloted aircraft are the static 
and dynamic characteristics that influence the ease and precision with 
which a pilot is able to perform the control task required in support of 
the aircraft mission flight phase. Thus, the handling qualities depend 
not only on aircraft characteristics and a mission flight phase but also 
on the pilot's subjective opinion of the ease with which he can perform 
the control task. 

To accurately assess the pilot's opinion of the handling qualities of 
an aircraft prior to first flight of a prototype, a groundbased simulation 
is usually required. In the early stages of the design, it is more eco- 
nomical to use a mathematical pilot modeling simulation because the design 
parameters can be easily adjusted. The pilot's assessment is then related 

to some scale such as the widely accepted Cooper-Harper pilot rating scale 
(Figure 1) . 

Much research has been done to determine the relations between the 
parameters of the rigid body, small perturbation equations of motion and 
the pilot rating. The handling qualities requirements for a rigid air- 
plane in Chalk et al. [1] are typical results of such research. Most of 
the airplanes in the past have been relatively rigid sucli that the elas- 
ticity of the airplanes do not contribute significantly to the pilot 
perceived handling qualities. 
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Recent advances in control -configured vehicles design and active con 
trol technology makes it possible to increase aircraft size and the utili 
zation of lighter structures in future designs. The elastic behavior of 
these vehicles is therefore becoming an appreciable influence in their 
handling qualities. Because of the potential adverse effects of mode in- 
teraction with the rigid body dynamics, there is a need for handling qual 
ities assessment in the preliminary design phase of new airplanes. 

General Background 

It is known that static aeroelastic deflections of an aircraft struc 
ture modify the aerodynamic pressure distributions which results in sta- 
bility derivative changes associated with the rigid body, small 
perturbation equations of motion. Early attempts to account for aero- 
elastic effects on aircraft stability and control took the approach of 
making static aeroelastic corrections to the aerodynamic stability de- 
rivatives [2-4] . The drawbacks of this approach as pointed out by Milne 
[4] are that in calculating modified stability derivatives one is to 
imagine the major parts of the airplane to be kinematically constrained 
at various points which do not have any real physical meaning, and if 
the overall -motion frequencies are of the same order as the lower 
typical vibration natural frequencies of the structure then the approach 
is invalid. 

For flying in high dynamic pressure environments, such as terrain 
following in turbulent air, the dynamic effects of flexibility are im- 
portant enough that they must be included as additional degrees of free- 
dom. A common approach has been to approximate the dynamics by a trun- 
cated set of superimposed orthogonal vibration modes. In this case 
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the phenomena of most interest are the effects of aerodynamic coupling 
between the various elastic modes and between elastic and rigid body 
modes, as well as elastic mode interaction with the feedback control 
system. Reference [5] was one of the earliest comprehensive studies 
of this problem. Hie most recent comprehensive work done under the 
AFFDL sponsorship is documented in reference [6] . 

The subject of handling qualities requirements and criteria for 
highly elastic airplanes in turbulent and high dynamic pressure environ- 
ments has been largely ignored. Much of the research on handling quali- 
ties has been concerned with relatively rigid, tactical military aircraft. 
Hie handling qualities parameters, such as phugoid, short -period, dutch- 
frequencies and damping ratios, which have been determined pertinent 
for such airplanes, are mostly meaningless for a flexible airplane with 
elastic mode frequencies close to the rigid body frequencies. When multi- 
ple frequencies are in proximity to one another, the pilot cannot easily 
discern individual modes of motion; rather his opinion of the transient 
dynamics will likely be based on the time history of the total motion. No 
performance criteria suitable for handling qualities specification are 
presently available for such higher-order responses. This is all too evi- 
dent in that no useful discussion of aeroelastic effects is included in 
the revision to the military aircraft handling qualities specification [1], 
It contains only the following statement: 

Since aeroelasticity, control equipment, and structural dynamics 
may exert an important influence on the airplane flying qualities, 
such effects should not be overlooked in calculations or analysis 
directed toward investigation of compliance with requirements of 
this specification (p. 497). 
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The specification is concerned only with desirable ranges of values 
on rigid body static and dynamic response parameters. There are methods 
available for estimating static acroelastic corrections to rigid body aero- 
dynamic stability derivatives [7]; however, the specification then re- 
quires the use of these in rigid aircraft equations of motion. It seems 
quite possible that the desirable ranges of parameter values could be 
significantly affected by elastic mode degrees of freedom, particularly 
when sane of the modes have natural frequencies of the same order of 
magnitudes as the frequencies of the rigid body alone. It is not at all 
clear that the handling qualities should be specified by rigid body dy- 
namic parameters when such mode interaction is present. In fact, the 
pilot could not tell, for example, how much of a given pitch angle re- 
sponse to command input is due to rigid body and how much to low fre- 
quency elastic modes. 

The key in developing handling qualities criteria and eventually 
specifications for severe mode interaction situations is to establish 
when and under what conditions the pilot can visually separate the rigid 
body response from the total response. In conditions when he cannot, a 
structural mode suppression control system probably will be required. 

Objectives and Scope of Study 

The primary objective was to develop an analytical method to de- 
termine the boundary between when the pilot can visually separate the 
rigid body motion from the total motion and when he cannot in terms of 
the small perturbation equations of motion parameters. This study is 
an extension of the experimental work done in reference [8], where a 
ground-based pilot- flown simulation was studied. Hie mathematical 
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pilot modeling simulation approach is used to assess the effects of 
mode interaction on the pilot opinion rating. 

An extension of the optimal control model for the human pilot [9] 
is made so that the effects of mode interaction can be assessed. TTie 
extension is motivated by an observation of the experimental evidences 
of reference [8] . 


Plan of Presentation 

A summary of past results is presented in Chapter II. The long- 
itudinal equations of motion for a flexible airplane are developed in 
Chapter III. The general description and flight condition of the flex- 
ible airplane under study are also described in that chapter. Numeri- 
cal values of stability derivatives for the equations of motion are 
given in Appendix A. The pilot modeling and its extension is presented 
in Chapter IV. Derivations of some singular perturbation techniques 
needed in Chapter IV are summarized in Appendix B. The major results 
are presented in Chapter V. TTie conclusions and recommendations appear 
in Chapter VI. 



CHAPTER II 


PAST RESULTS 

The only research of which we are aware that is directly relevant to 
the subject is documented in Crother [10] and Yen [8] . 

The results of Morth American Rockwell in Crother [10] were for an 
early version for the B-l aircraft and included piloted simulator evalu- 
ations of tracking performance in turbulence. They concluded that the 
structural dynamics appeared as essentially a nuisance oscillation to the 
pilots and did not significantly effect tracking performance. However, 
the longitudinal dynamics of their configuration were very close to Case 
1 of our results. Thus, it is not surprising that the elasticity did not 
significantly degrade pilot opinion; it was merely a ripple on the rigid 
body response. 

In the work of Yen [8], the effects of parametric lowering of the 
undamped natural frequencies of the first two symmetric elastic modes 
of a flexible aircraft were investigated. A pitch tracking task, which 
included phugoid and short period dynamics, was programmed on a fixed- 
base simulator with a CRT attitude-director display of pitch command, 
total pitch angle and pitch error. The display and its variables are 
depicted in Figures 2 and 3. 

The attitude -director equations are 

= 0(t) ' O- 025 ^^) - 0.029s 2 (t) (3.1) 

pitdi error = e 0 = (3.2) 
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where 0,025 and 0,029 are the slopes of the two elastic modes at the pilot 
station. The flight condition is Mach 0,85 at sea level density, 

Four pilots each flew eight cases which were combinations of clastic 
mode interaction. The cases arc shown in Table I, Case 1 is the original 
dynamics. Case 6 was the most difficult, where the free-free elastic mode 
frequencies were set at 6.93 rad/s for both modes. This resulted in the 
phugoid mode splitting into positive and negative real roots. Sample time 
histories of one pilot's tracking difficulty on Case 6 are shown in Fig- 
ure 4. Note the large amplitudes of the elastic modes' contributions to 
total pitch 0 i relative to rigid pitch e. This made it very difficult 
for the pilot to visually separate rigid from elastic pitch. The average 
of the four pilots' ratings of Case 6 was 6.7 on the Cooper-Harper scale. 
Contrast this with a 1.6 rating on Casel, the original dynamics. This 
work has clearly established the potential seriousness of elastic-rigid 
body low frequency mode interaction to handling qualities specifications 
and pilot rating. 
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Figure 3. The Airplane Attitude Corresponding To Above EADI 



TABLE I 


NATURAL FREQUENCIES AND DAMPING RATIOS OF EIGHT CASES 


Case 

# 

“l 

rad/sec 

L 

2 S sp 

rad/sec 

“sp ^ph 

rad/ sec 

> 

rad/sec 

? le 

w le 

rad/ sec 

C 2e 

w 2e 

rad/sec 

P.Q.R. 

1 

13.59 

21.18 

0.5339 

2.806 0.0197 

0.0708 

0.0494 

13.312 

0.0215 

21.354 

1.6 

2 

9.17 

21.18 

0.5235 

2.5724 -0.00060267 0.0573 

0.08769 

8.7891 

0.0213 

21.356 

2.0 

3 

6.16 

21.18 

0.5217 

1.7691 

Real Roots 

+0.090978 

-0.076723 

0.1999 

5.8669 

0.0213 

21.357 

5.9 

4 

13.59 

4.79 

0.6872 

1.5745 

Real Roots 

+0.14654 

-0.13167 

0.05284 

13.270 

0.1137 

5.9702 

3.1 

5 

11.66 

11.66 

0.5436 

2.5819 -0.0001122 0.0537 

0.0773 

11.801 

0.0162 

11.574 

2.0 

6 

6.93 

6.93 

0.7028 

1.3665 

Real Roots 

+0.17581 

-0.15307 

0.1919 

7.3305 

0.007599 

6.9178 

6.7 

7 

10.25 

9.75 

0.5517 

-0.0483 

0.0282 

0.1129 

10.234 - 

•0.0004277 

9.8978 

2.3 

8 

10.68 

9.27 

2.3893 

-0.0541 

0.0256 

0.11021 

10.347 

0.0005306 

9.7781 

1.9 





CHAPTER III 


EQUATIONS OF MOTION 

The equations of a flexible airplane consist of an overall spatial 
motion, that is, a rigid body motion, and a local deformation due to 
its inherent flexibility. The basic principles underlying the equations 
of motion are the conservation of the linear and angular momenta and an 
internal equilibrium due to elastic deformation. The equations are 
written in terms of the x,y,z body-fixed axes frame of reference. The 
orthogonal axes are chosen such that the x-axis passes through the cen- 
ter of gravity of the airplane and points forward parallel to the free 
stream steady-state (trim) velocity, the y-axis points out to the right 
wing, and z-axis points downward. 

The deformation of the airplane is expressed in terms of natural 
mode shapes and generalized coordinates. The airplane is assumed to 
be a plate like structure in the normal mode (in-vacuum vibration modes) 
calculation. 

In the next four sections the pertinent assumptions and the equa- 
tions of motion are summarized. The detailed derivation and related 
discussions can be found in [ 11 - 14 ]. 

Small Perturbation Equations of Motion 

For a cruise, level flight condition at a trim speed U , the small 

o 

perturbation longitudinal equations of motion are given by: 
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mu = -mg9 + X u (u + u ) + XJa + a ) + X 6 <5 

b ^ c 

m u 0 (<S - 9) = Z u (u ♦ u„) ♦ z a (a * a s ) + z^d . y ♦ Zj( j + y 

* Z. 9e +.?, [H_f. + 12_> .1 

(Se 1=1 L 3£ ^1 3^i lJ 

I y e = M u (u + u g ) + M a ( a + a ) + M.(4 - q ) + M- (0 + q ) 


(3.1) 


+ M. fie + .E, [|^e. + 

6e i=l l aTpi ^|7V 


3M 


, 3Q «- 

”i [ <i * 2! i“i«i * “i ? i) " C ~5£-^ ) (« * « g ) * -Sqi (9 * q g ) 


where : 




3Q f 


>, <?• 3Q f. 

* C M7 )S e * ^<4 - V * A ^ «k 'jT'k 1 ' i=1 ' 2 - 


s «k 


n s d( )/dt 


u Perturbation forward speed 

u g Perturbation of u due to gust 

0 Perturbation pitch angle (rad) 

q Perturbation pitch angle rate (rad/s), (q=t>) 

q g Perturbation of q due to gust 

a Perturbation of angle of attack (rad) 

a Perturbation of a due to gust 

o 

Perturbation of elevator deflection (rad) 
q Generalized coordinate of the i th elastic mode 
m i Generalized mass of the i th elastic mode 
Generalized force of the i th elastic mode 
u i in-vacuum elastic mode undamped natural frequency of i*^ 1 mode 
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The equations (3.1) are rewritten in a state-variable formulation 
as: 

*a (t) * Va (t) + B a u a (t) + E a w a (t) C 3 - 2 ) 

where 

x a (t) = col. [u g , o , a g , q g , u, a, 0, 0, C 2 ,..., 

U ft) = 6 
a v J e 

Hie state a and the input vector w (t) are to be discussed in the 
fal a 

next section on the turbulence model. 


Turbulence Model 

The turbulence model is derived from the Pryden gust power spectra 
which have the forms [1], 


♦u « ■ °u — T— , 

S o 1 + f \x ^ 2 




l 1 * 3( ir 

, c ^ 2 w o 

V B) " w ^ ^ — r. 


<f> («) = ~T~ Vo) 

« U o 8 


HJ J 


4> a (oj) = 
4 g 


i + ( 4 !iS -) 2 

V 7T U 


+ w C«) 


where 


2 1 


2? 4 do, i = U g , w g , cx g , and q 
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a temporal frequency (rad/s) 

U Q true air speed 
b w wing span 

L u , 1^ gust scale factors which depend on the altitude 

The time domain representation of the turbulence as a shaping filter 
with zero mean, gaussian, white noise processes and ^ as inputs is: 
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and 



Attitude Director Equations 

The total pitch angle time history that the pilot feels and sees, 

either on the outside horizontal or the attitude indicator display, is 
given by [10] : 


y a^ x p ,t:) ' 0( ^ t ) ‘.2 KOO £,(t) 

i=l J 1 J 


(3.5) 


= 0(t) - e e (Xp,t) 

where x indicates pilot fuselage station, <K(xJ the slope of the 
• th J P 

j symmetric elastic mode at that station, e(t) the rigid body pitch 

angle, and 0 e (t) the elastic contribution to the total pitch angle 
(Figure 5) . 

The equation (3.5) can be written in terms of the state variables 
defined in equation (3.2) as follows: 


y a^ » C a x a (t) 


where 


C a = col. [0,0,0, 0,0,0,1, 0,-*J-*' 0,0,.. .J 



;m=» 
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B-l Flight Condition 

The B-l bomber was chosen for this study because it exemplifies the 
trend toward more elastic structures for future large aircraft f* The 
total length of the B-l is 46 m (151 ft.). The reference wing span util- 
ized at the flight condition in Table II is 41.7 m (136.67 ft.). The 
values of the stability derivatives and the necessary data for the equa- 
tions of motion are given in Appendix A. 


TABLE II 

B-l FLIGHT CONDITION 


Mass = 103,370.15 kg (7085.0 slugs) 

Mach No. =0.85 

Velocity = 289.4 m/s (949.0 fps) 

eg at fuselage station = 40.67 m (1061.2 in) 

I y = 8.0 x 10 6 kg-m 2 (5.9 x 10 6 slug - ft 2 ) 

S w = 180.8 m 2 (1946.0 ft 2 ) 

c w = 4.67 m (15.33 ft) 

b w » 41.7 m (136.67 ft) 



CHAPTER IV 


PILOT MODELING 

The human pilot in a manual control task can be modeled as an 
active feedback element in the aircraft control system. The quasi- 
linear model and the optimal control model (001) are the two models 
widely used in this way. The quasi- linear model has the analytical 
description in terms of the frequency-domain control system design 
technique, while the optimal control model is based on the time- 
domain or optimal control theory. Since the analysis in this study 
is mostly in the time-domain, the optimal control model is employed 
through out. There are other reasons for employing the optimal con- 
trol model which will be discussed later. 

The optimal control model of the human pilot was originally 
developed by Kleinman, Baron and Levison [9], The fundamental assump- 
tion underlying the OCM is that the well -motivated, well-trained human 
pilot will act in a near optimal manner subject to the pilot's intern- 
al limitations and understanding of the task. By specifying human 
limitations, the optimality assumption gives a model that adapts to 
task specifications and requirements automatically and not through 
a subsidiary set of adjustment rules as has been done in the quasi- 
linear model. Thus, for a new situation, the optimal control model 
can be modified by just determining the operative limitations and the 
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new control task. The review of the past applications of the model 
can be found in reference [15]. 
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Model Description 

The structure of the optimal control model of human pilot response 
is shown in Figure 6. The aircraft dynamics, which also include noise 
shaping filters of the turbulence, are described by the linear, time 
invariant equations. 

*a (t > = Va™ + * W'J (4.1) 

x a (0) = given 

where 

x a (t) = aircraft and shaping filters state vector of dimension Na 
u a (t) = pilot's control input vector of dimension Nu 
w a (t) = disturbance vector of dimension Nw, each of which is an 
independent zero mean, Guassian white noise process with 
covariance 

E {w a. (t) w a.^ } = W i 6(t ' a) ’ 1 = 1.2,..., ^ (4.2) 

The display variables are given by a linear combination of state 
variables . 

y a Ct) = C a x a (t) (4-3) 

where 

y a (t) = displayed vector of dimension Ny 




Figure 6. Optimal Control Model of Human Pilot Response 
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The usual assumption in the model is that if a quantity y^ is explicitly 
displayed to the pilot, he also derives the rate of change y^. Thus, 
y a (t) contains both position and velocity information of a displayed 
signal, but no higher derivative information. 

Internal Model 

The information processor in the pilot model operates on a noisy, 
delayed version of the displayed variables to obtain a "best" estimate 
of the aircraft state vector. This is accomplished by a Kalman filter 
and a least-mean-square predictor and makes use of an internal model. 

The internal model of the pilot may be considered to consist of [16] . 

1. Knowledge about the overall behavior of the aircraft 
under control and about the possibilities to control it. 

2. Knowledge about the disturbances acting on the aircraft 
and the way they will influence it. This knowledge 
will be of a statistical nature. 

3. Knowledge about the task to be performed. 

In many instances , the assumption that the internal model is an exact 
replica of the system model, i.e., perfect internal model, appears to 
be a satisfactory one [17] . There are situations in which the assump- 
tion of a perfect internal model does not appear tenable. In a highly 
complex system, i.e., one with a large number of state variables, with 
a single display it is unlikely to be modeled perfectly by the pilot. 

It is of interest in this study to determine the effects of high 
frequency oscillation, contributed mainly by the elastic modes, on 
the handling qualities and pilot rating. From a past experiment with 
ground based simulation of the elastic airplane [8] , the pilot action 
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in a pitch tracking task closely resembles that of the rigid body 
pitch error rather than the total error displayed to him. This shows 
the pilot's ability to filter out the high frequency oscillation in the 
total pitch response and leads to a hypothesis that the pilot uses the 
slowly varying dynamics subsystem as the internal model. 

Hie decomposition of the slowly varying dynamics from the aircraft/ 
disturbance dynamics, Eqns. (4.1) and (4.3), is accomplished by the 
singular perturbation technique (Appendix B) . The aircraft/disturbance 
dynamics can be written in the form: 

*l (t) = + A i2 x 2 (t) + B l u a (t) + H l w a (t) (4.4) 

yx 2 (t) = A 21 Xj(t) + A 22 x 2 (t) + B 2 u a (t) + l : . 2 w a (t) 

y a00 = C 1 x 1 (t) + C 2 x 2 (t) 

where 


~ rigid body and noise shaping filters state vector 

x 2 (t) = elastic modes state vector 

y = a small positive parameter which arises due to the presence 
of high frequency elastic modes and can be an unknown in 
this analysis. 

Then, by letting y->0 + we get 

*d« = A d V 1 ) + B d u a Ct) * E d w d'« C4 S) 

y d (t) ” C d x d (t) * D d u a (t) 

where 
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By using an imperfect internal model, the computational task be- 
comes formidable since it involves the solution of the matrix delay 
differential equation [17] . The structure of the original optimal con- 
trol model can be modified from Figure 6 to that of Figure 7. This mod- 
ification does not affect the prediction capability of the model very- 
much as has been shown in references [17-19]. In this modified model, 
the pure time delay x is approximated by a first-order Pade polynomial. 
Thus, the relation between i^(t) and r(t) in Figure 7, which was origi- 
nally expressed by 


u a Ct) = r(t-x) 

or in the Laplace transform operator s, 
u a (s) = e' TS r(s) 


is now approximated by 


V*) 


~ ~ s + 2/x. 


s + 2/x 


which can be expressed in the state variable form as 


(4.6) 

(4.7) 

(4.8) 


u a (t) = z(t) - r(t) 
where 


(4.9) 


z(t) =-|- I z(t) + ^ I r(t) (4.10) 

I is an identity matrix of appropriate dimension. The time delay x 
is normally 0.1 to 0.2 sec. 


Human Limitations 


Other than the time delay, the pilot has inherent limitations of 
perceptual noise and perceptual indifference thresholds on displayed 
information. The time delay has been compensated for in the control 
action as shown in Eqns. (4.9) and (4.10). The other quantities are 
associated with the observation process in the pilot model, so the 



Disturbance 




Figure 7. Modified Optimal Control Model of Human Pilot Response 
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pilot is assumed to perceived v ft) • 

1 V a V t} whlch 15 a noisy version of y (t), 
1 ^ 0 • ^ 3 


V tJ = y a ft ) + yo 


(4.11) 


where the observation threshold is replaced by the Random Input Describing 
unction N(o) and incorporated in the observation noise, v (t, is . zero . 
*an, gaussian, white noise process with autocovariance. 

t* r_ r » 




When directly viewing y a (t) , the associated COTariance 

1 w y 

7 i 


(4.12) 


>i 


e2 

a. 


(4.13) 


where 


-1 


°i " a yi N(a y i ) 

a y. (t)f 

i a i 

"V = erC '(~g, describing function gain of threshold 

erfc = error function 

a i = half widt h of dead zone element 

- noise/signal ratio at full attention on indicator i 

- O.Oli or -20 dB normalized power density level 
fj - attention allocation to display indicator i 

For the total of X indicators, neglecting the time spent in interinstru- 
nicnt scanning, we have 


i=l f i = 1> 0<f.<l 

value of 
is minimized. 


(4.14) 


isUntod/ 1 15 Ch ° Sen ^ that ^ C ° St functioral of 'he pilot model 
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Task Definition 

The important assumption about the optimal control pilot model is 
that the pilot's control task is adequately reflected in the choice of a 
control r(-) that minimizes the cost functional of the form 

J(r) = li»E{ | /J (y'CtXyjt) + r'(t)dj ( r(t)]dt) (4.15 

conditioned on the perceived information y p (-J. is a specified con- 
stant, symmetric, nonnegative definite matrix which depends on the task 
specification. The control rate term is used to account for the pilot's 

Imitation on the rate of control motion and introduces first-order neuro 
muscular dynamics in the pilot model. 

Hie selection of the weighting 0 = diag. [q ] is such that 

7 y i 

q - 1 2 

y i y pi, max (4.16) 

where 

y Pi,max 1S the maximum desired or allowable value of y Unlike Q , 
the weighting Q R - diag. [q r ^], a positive definite matrix, is not speci- 
fied before the pilot model equations are solved. It can be shown that 
the pilot control law which minimizes (4.15) takes the following form: 
T n r(t) = -r(t) * V t) * v B (t) (4.17) 

Hie matrix T n is assumed to be in the fonn T n = diag. [t n ] , i=l,2,..., Nil. 
Hie scalars are a neuromuscular time constant of humai limbs, which 
has a typical value of 0.1 sec., independent of the system to be control- 
led. Thus, the weighting q^ are adjusted iteratively until each 0.1 
sec. If the resulting q^ weighting is such that 1/^ is much 1 
greater than the physical rate at which one can move 1 
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control r i? then q r = |l/f\ (t), max | 2 must be used. Though, this 
rarely happens- except for highly unstable aircraft dynamics or an 
aircraft flying through very severe turbulence. 

The motor noise v m (t) is a zero-mean, gaussian, white noise 
process, with autocovariance 

Etv m w V 0 ))- V (t -°) 

and V^is known to scale with E(m 2 (t) } ,i.e., 

where the typical value of the motor noise/signal ratio 
?m.=0.003ir. 

Equations (4.9) and (4.10) may now he augmented to Eqn. 

(4.1) to define an augmented system of equations 

*c (t: > = A c x c (t) + B c u m (t) + E c w c (t) 

: C X C ( 


(4.18) 

(4.19) 


(4.20) 


y a (t) = Cx r (t) 


where 


x c = f x a * z * r l 


\ B -B 

3 3 3 


A = 
c 


0 -2/t 4/t 
ID 0 




B c =[0, 0. T/ 1 ] 
C c -[C a ,0,0] 


re a o 1 


E = 
c 


0 0 


0 T 


-i 


n - 1 


w c “I"' V 


Equation (4.20) is the "actual" dynamics to be controlled by the pilot, 



30 


The Pilot Model 


The pilot's control input r(t) that minimizes J(r) is generated 

based on the augmented system of the internal model (4.5) and the delay 

compensation (4.9) and (4.10), i.e., 

x (t) = A x (t) + B r(t) + E w(t) 

s (4.21) 

y s (t) ■ C s + D s r(t) 

where 



4/t ^ 



[C d , Dj] 






The command control of the pilot is given by 


U n> (t) ■ - L opt V 1 ’ 
- -L*x t (t) 


(4.22) 


L * * l L opt 01 

■ M 

Lx J 

L opt = P 22 Qr 


where 
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T n • P 22 P 12 


P = 


P 11 P 12 


L P 12 P 22J 


satisfies the equation 


A ; p * PA o ,c ; < )yV PB o ( !R 1 »; i, = » 


(4.23) 


where 


A o * 


A s B s 


0 0 
B o - $ 

C o ■ I c s D sl 

Hie state x £ (t) is the best estimate of x t (t) generated by a Kalman filter 
*t ( t) - + * KtXp - C 0 i t (t)] (4.24) 

where 


K = EC'V 
o y 


(4.25) 


and z satisfies the equation 


V + 


Ik' + EWE' + ZC'V 
t o t o o y 


-1 


C o E 


(4.26) 


where 
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V 


IV 0 

o v_ 


m 


Combining equations (4.11) with (4.20), (4.22) and (4.24) yields the 
closed loop system 


= • B C L * n c “ c Ct) 

x t (t) = (Aj - BjL)^ * K[C c x c - C o x t ♦ V (t)] 


or 


= Fi|» + Gw 


where 


<l> = 


F = 


L KC c 


-B C L* 
A, 


-BjL* 


-KC 


G = 


w = 


w 


L v yj 


Thus, 


cov V> = 


cov x x' 
c c 

cov x.x' 
t c 


COV XX* 

c t 

COV X t X' 


= V 


is the solution of 

V = FY + W'+ GOG' 

where 


o - 


W c 0 

0 V 


(4.27) 


(4.28) 


(4.29) 


(4.30) 
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Hess [18] has formulated a pilot 

™ pu<>t „ dellng Procedur ; ztz tcc,miquc for - “ 

validated in a variety Qf ^ ^ ' J“ haS ^cesstuUy 

stated as follows: 6 ratinfi technique can be 


(1) the indes of perforce and model 
optimal control pilot mnH i • rajneters in the 

e ing procedure yield a dvn 

sentative model of the human pilot, lly repre- 

(2) the variables selert^ r 

performance are directly h 0r lnClUSl ° n ln 016 index of 

directly observable by the pilot, 

(3) the weighting coefficients • 

™ *«n as the squares of the ‘ " ° f 

deviations of the respect' reci P roc als of maximum "allowable" 

sonent withdre ~ *- -iations are con- 

task as perceived by the pilot. 

Then 

the numerical value of the indr>v r 

the modeling procedure can be related to a Perf0I * nCe rcSUlti " 8 f ™ 
- Pi- assigns to the vehicie I tl C “ — — 

n a n ' 


where 


P ’ 0,R * “ 2 - 51 ^ (10 J) + o.3 


(4.31) 


P j’ R ’ _ = Pil0t ° Pini0n ^oper-Harper scale 

value of the perfonnance index 

Computational Algorithms 

There are two minr 

rne response variables 
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for a piloted-aircraft manual control task. A digital computer program 
STDOCM is a modification and extension of the program PIREP [21] written 
for operation on CDC-6600 at Wright Patterson Air Force Base to implement 
the standard optimal control model of the human pilot. A program MODOOI 
is developed to implement the modified optimal control model of the human 
pilot which is presented in this chapter. Both programs are written in 
Fortran IV for operation on IBM system 370/168 at Oklahoma State Univer- 
sity and are available from Professor R. L. Swaim, School of Mechanical 
and Aerospace Engineering, Oklahoma State Unversity. 



CHAPTER V 


EFFECTS OF ELASTIC MODES INTERACTION 
ON HANDLING QUALITIES 

The elastic modes interaction with the rigid body dynamics is intro- 
duced to the large flexible aircraft model by parametric lowering of the 
undamped natural frequencies of the two elastic modes. This will cause 
controlling the rigid pitch angle by observing the total pitch error to 
become more difficult as indicated in the past experimental results [8] . 
The standard optimal control model for the human pilot and modified model 
presented in Chapter IV are applied to the illustrated cases of varying 
elastic modes interaction. It will be shown that the standard optimal 
control model gives the misleading results when there is a severe modes 
interaction between the elastic modes and the rigid body dynamics. The 
modified model gives more consistent results with the experimental data 
on the effects of elastic modes and the rigid body dynamics. The modified 
model gives more consistent results with the experimental data on the ef- 
fects of elastic modes interaction on handling qualities and pilot ratings 
than that of the standard optimal control model. In this chapter the il- 
lustrated cases used in the computer simulation study are described. Then 
the simulation results are presented. Finally, the separation boundary 
which can be used as an indicator of when the pilot can or cannot visually 
separate the rigid body motion from the total motion is presented. 
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The Illustrated Cases 

The ten illustrated cases are obtained from the equations of motion 
(3.1) in which the natural frequencies of the two elastic modes are para- 
metrically reduced. Dynamic characteristics of each case are specified by 
four modes: phugoid mode, short period mode, elastic mode 1, and elastic 

mode 2 as shown in Table III. These ten cases will exemplify most of the 
situations in which the handling qualities and pilot ratings would be af- 
fected differently by the two elastic modes included in the model. 

The lowering of the elastic mode natural frequencies resulted in mode 
interaction which lowered the coupled short-period and phugoid frequencies 
or made one of them split into positive and negative real roots. A full 
state feedback control law is used to place the roots of the characteris- 
tic equation at precise values for each case. The rigid body dynamics are 
maintained to be the same as Case 1 and the elastic mode coupled frequen- 
cies placed at original values before the state feedback control law was 
applied. This will ensure that the pilot ratings are based on the rela- 
tive amplitudes of rigid and elastic pitch angle responses and not on poor 
rigid body dynamics. 

The simulation results on the ten cases by using the standard 
optimal control model (OCM) for the human pilot and the modified model 
are shown in Tables IV and V, respectively. These results clearly indi- 
cate that when there are severe modes interaction, such as Cases #3, 6 
and 9, the standard OCM gave very low pilot ratings predictions which are 
inconsistent withthe experimental results [8], In contrast the modified 
OCM gave more consistent results since it includes the visual separation 
process of the rigid body response from the elastic modes response which 
the pilot has to accomplish when the amplitude of the high frequency 



TABLE III 


NATURAL FREQUENCIES AND DAMPING RATIOS OF TEN CASES 


Case 

If 

rad/s 

u)2 

rad/s 

°ph 

rad/s 

Sh 

—T$ 

rtS/s S P 

-le 

rad/s 

5 le 

w le, 

rad/s 

2e 

1 

13.59 

21.18 

.0665 

.0312 

2.9334 .5209 

13.236 

.0497 

21.395 

.02112 

2- 

8 

21.18 

.04614 

.001376 

2.581 .4992 

Real Roots 

7.508 

.1127 

21.390 

.02104 

3 

4 

21.18 

.1412 

.1336 

+ 1.652 
- 2.266 

4.544 

.3892 

21.380 

.02102 

4 

13, 59 

15.00 

.06345 

.02899 

2.889 .5247 

13.04 

.04508 

15.480 

.03052 

5 

r 

8 

15.00 

.04489 

.001946 

2.586 .5031 

Real Roots 

7.400 

.1073 

15.340 

.02787 

0 

4 

15../0 

.1411 

.134 

+ 1.568 
- 2.172 

4.356 

.3990 

15.320 

.02761 

7 

13.59 

13.59 

.06203 

.02797 

2.87 .5262 

12.76 

.03821 

14.380 

.02923 

8 

8 

13.59 

.04433 

.002208 

2.588 .5048 

Real Roots 

7.345 

.1042 

13.990 

.03117 

9 

4 

13.59 

.1411 

.1350 

+ 1.527 
- 2.126 

4.269 

.4035 

13.970 

.05062 

. 10 

8 

8 

.03801 

.005574 

2.608 .5245 

6.403 

.04999 

9.390 

.08138 



TABLE IV 


PERFORMANCE PREDICTION BY THE STANDARD OCM 


CASE 

# 

0 • 
rms 

(deg) . 

• 

0 

rms 

(deg/s) 

6 

rms , 
(rad x lO'j 

P. 0. R. 

1 

.3895 

1.650 

1.143 

1.0 

2 

.5422 

2.343 

2.475 

2.6 

3 

.5042 

1.573 

2.234 

1.0 

4 

.3612 

0.917 

0.910 

1.0 

5 

.5096 

2.335 

2.266 

2.5 

6 

.4604 

1.721 

2.061 

1.2 

7 

.3437 

1.127 

1.392 

1.0 

8 

.4955 

2.288 

2.142 

2.4 

9 

.4445 

1.750 

2.000 

1.2 

10 

.3833 

1.591 

2.810 

1.0 
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TABLE V 


PERFORMANCE PREDICTION BY THE MODIFIED OCM 


CASE 

§ 

9 rms 
(deg) 

6 

rms 

(deg/s) 

6 

rms , 
(rad x 10 5 ) 

P.O.R. 

1 

.3663 

1.848 

0.829 

1.2 

2 

.5724 

2.832 

1.199 

3.3 

3 

.8123 

3.721 

1.014 

4.7 

4 

.3157 

0.999 

1.024 

1.0 

5 

.5345 

2.746 

1.195 

3.2 

6 

.7374 

3.470 

0.872 

4.3 

7 

.3228 

1.428 

0.869 

1.0 

8 

.5172 

2.657 

1.181 

3.0 

9 

.7043 

3.360 

0.833 

4.1 

10 

.4187 

2.091 

1.167 

1.8 
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elastic modes is getting larger. On the other hand, when the modes inter- 
action are small, both pilot modeling techniques gave almost the same pre- 
dictions. 

The Separation Boundary 

The separation boundary is defined as the limit of when the pilot can 
or cannot visually separate the rigid body motion from the total motion. 

The visual separation is essential in controlling the rigid pitch angle or 
other rigid body parameter of the elastic airplane when only the total 
pitch angle response or a corresponding parameter is available to the pilot. 

From Chapter IV, two kinds of pilot modeling techniques are discussed. 
The standard OCM is the pilot model that asssumed that the pilot has the 
perfect internal model of the aircraft/disturbance dynamics. This model 
will give the best possible pilot opinion rating (P.O.R.) in any tracking 
task. The other model, the modified OCM, is the one with a slowly varying 
internal model. In this model, the pilot is assumed to be able to com- 
pletely separate the slowly varying or the rigid body motion from the to- 
tal motion. The difference between the P.O.R. 's, i.e., AP.O.R., of 2 is 
chosen to be a separation boundary. That is if AP.O.R. is greater than 
or equal to 2, the pilot cannot visually separate the rigid body motion 
from the elastic motion in the display. 

It is known that the pilot opinion rating depends on many factors 
such as the intensity of turbulence and the level of difficulty of the 
task. To study the modes interaction effect all other effects should 
be kept at their nominal values. That is without severe modes interaction 
effect the other parameters should be set such that P.O.R. is equal to 1. 
Once the P.O.R. has been initialized for some specific task, the separation 
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boundary can be found from the P,O.R. prediction of the modified 0Q4 
alone. This is becduse the standard OCM will give almost unity P,0,R. 
in the severe modes interaction cases, if the rigid body dynamics have 
been maintained at the known good handling qualities specifications. So, 
instead of using AP.O.R. = 2 as the separation boundary, the modified 
OCM's P.O.R. of 3 can be equally well used as the separation boundary 
provided the proper initialization mentioned above has been done. 



CHAPTER VI 


CONCLUSIONS AND RECOMMENDATIONS 

A model for the human pilot in a manual control task using the 
optimal control techniques has been developed for predicting the 
handling qualities or the pilot ratings of a large flexible aircraft 
where the elastic modes interaction with the rigid body dynamics is 
significant. The separation boundary concept, which will tell when 
the pilot can or cannot visually separate the rigid body motion from 
the elastic one, has been introduced. The techniques developed here 
make it easier to investigate the modes interaction effect on the 
handling qualities in a preliminary design stage before the first 
prototype has been built. If the handling qualities are severely 
affected by the elastic modes interaction with the rigid body dynamics, 
the elastic modes suppression control system should be designed and 
implemented along with the stability augmentation system of the rigid 
body dynamics. 

A comparison of the model predictions with the past experimental 
data shows that the modified optimal control model for the human pilot 
developed here is much better in predicting the elastic modes effect 
on the handling qualities than the standard optimal control model. 

This is due to the fact that the mode decomposition mechanism has 
been incorporated into the modified model. 
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However, only one longitudinal trim flight condition was investi- 
gated, its validity for different flight conditions can only be confirmed 
by conducting more investigations. 

Future work should also include the computational aspect of the 
model. In the pilot model parameters identification, it is required to 
solve an (n a + n s + 4)*(n a + n g + 4) matrix equation (4.30) for the mod- 
ified model instead of solving an (n + l)*(n„ + 1) matrix equation in 
the standard OCM, where n a , dimension of an aircraft/disturbance dynamics, 
is 12; n s , dimension of a slowly varying part of the aircraft/disturbance 
dynamics, is 8. This will risk the numerical instability when one has 
to include more elastic modes in the aircraft/disturbance dynamics. An 
alternative structure of the pilot model should be explored to ease the 
computational burden of the high order aircraft/disturbance system. 
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APPENDIX A 


NUMERICAL VALUES OF STABILITY DERIVATIVES 
AND EQUATIONS OF MOTION 


The low level penetration flight condition for the B-l was supplied 

by the B-l System Program Office at Wright -Patterson AFB from Rockwell 

International unclassified documents, The stability derivatives used in 

euqations (3,1) were based on preliminary aerodynamic analyses, but closely 

representative of the vehicle that has been flying. The relations between ' 

dimensional force and moment and elastic force derivatives as a function 

of nondimens ional stability derivatives are given in Table VI and VII, 

respectively. The non-dimensional stability derivative values for the 

unqugmented vehicle are given in Table VIII. The gust specifications for 

the study vehicle are given in Table IX. Finally, the A^, A-^i 

B , B 2 /H, Cp C 2 and E^ matrices for the unaugmented airplane are 

given in Tables X to XVIII, respectively. The matrices F^ and D are 

zero matrices and the matrices A , B , C and E are defined as follows: 

a ’ a ’ a a 


A a s 


A 11 A 12 


A 21/p A 22/u 


B, 


B a S 


B 


2/y 


C = [C. C 7 ] 
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TABLE VI 


nrmsioNAL force and moment derivatives 

AS A FUNCTION OF NON-DIMENSIONAL 
STABILITY DERIVATIVES 


a x = plCS/2 


a~= a.c/2U 
2 1 o 


X = a.c /U 
u lx' o 
u 


Z u - a l C z /l 'o 
u 


M . = a lC C /U 
u I m ' o 


X = a-C 
a lx 

a 

Z = a.C 

a 1 Z 

a 

M = 

a 

a.cC 
1 m 

a 

X. = a 7 C 
a 2 x. 

a 

Z. = a-C 

01 Z Z • 

a 

M. = 
a 

a 7 cC 
z m. 

a 

X 8 ’ a 2 C x. 

6 

Z ‘ ■ a 2% 

M 0 " 

a 7 cC 
2 m 6 

X, - a.C 
S e 1 

V a ^6 c 

M 6 = 
e 

a.cC 

X = 1 
6 e 

Z* = a,C, 

<5 1 Z, 

e 6 e 

M « = 
e 

a i cC 
1 m 6 e 

X, = 1 

6 t 

Z, = 0 
s t 

M r = 

6 t 

0 

V a l\ 

Z £ = a l < 'z 
1 E l 

V 

1 

a.cC 
l m 

h - a i V u o 

i ^ 

V a iV u o 
1 1 

M. - 
1 

a.cC /U 
1 m-' o 

1 

X 5 = a l C x 

2 

Z C = a l C z 

2 S 

M 5 “ 
2 

a-.cC 

1 

2 

X 5 ’ a l C x/ U o 
2 5 

2 

Z « ■ h C z/ V o 

2 S 

M • = 
S 

2 

a n cC /U 
1 o 

2 
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TABLE VII 


DIMENSIONAL ELASTIC FORCE DERIVATIVES 
AS A FUNCTION OF NON-DIMENSIONAL 
STABILITY DERIVATIVES 


a l " p *o a 2 = a i c /^ LI C 


Q 5 a i c c 


Q 5 = a 1 C^ 

o o 

-a 


a 2 C ? 


Q ? - a 2 C ? 

2 d 2 d 


<k =a 2 C ? 


% , = a 2 C £ 


" a^ 

\ ‘S 


Q c - a jC^ 

2 e 2 c 


Q C . = a l C C / U o 
2 2 


Q C . ■ a l C « /"o 

2 C 2 C 

2 2 


2 2 


Q c - ajC^ 


*5 . = a l C 5 / U o 


l S 16 


Q C = fl l C 5 
2 6 2 6 
e e 
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tabu: viii 


STABILITY l)i:iUVATIVi:S BOR B-I BGMBIR 
IN MAGI 0.85 BLIGHT CONDITION 


c x = 

-0.08066 

c 2 - 

-1.9659 

C = -0.4546 
m u 

x u 


Z u 


C = 

-0.08500 

C = 

-3.9367 

C = -1.41052 

X 

Z 


m 

a 


a 


a 

c = 

0 

c = 

-5.0 

C = -11.005 

X. 

Z. 


m. 

a 


a 


a 

c = 

0 

c = 

17.8558 

C = -35.7556 


z e 



c = 

0 

c = 

-0.9426 

C = -2.799 

X 6 

Z 6 


m 6 

e 


e 


e 

C X = 

0 

C z = 

-0.02922 

C = -0.0348 



5 


5 

1 


1 


1 

0. = 

0 

c 2 . = 

-0.6592 

C = -1.32169 
m. 

x 4 


z i 


i 


i 


1 

C x = 

0 

c, = 

0.015 

C = 0.03787 

K 


2 ? 


2 


2 


2 

C v = 

X 5 

0 

V 

0.4733 

= 1.233 

2 


2 


2 

C C = 

-0.06478 


0.48975 




'la 




0.02469 

n 

II 

0.48779 




id 



C C = 

-1.47658 

» 

3.97547 




% 



C 5 = 

0.00064 

c « ■ 
2 t 

0.00451 






i 


1 





51 


= -0.07243 

H 

1 

C ? = -0.0014 

2 

C = 0.0765 

^ . i 


= -0.19635 


TABLE VIII (Continued) 
= -0.07333 

n 

1 

= -0.0051 
2 C 

2 

C = -0.2588 
2 

= 0.3939 


TABLE IX 

GUST SPECIFICATIONS 


Parameters 


"Value 


w 


u 


w 


u 


U r 


'w 


6 fps 
10.8 fps 
300 ft 
970 ft 
979 fps 
136.68 ft 



TABLE X 


A n MATRIX 


-.9777 

0 

0 

0 

0 

0 

0 

0 

0 

-3.1633 

0 

0 

0 

0 

0 

0 

0 

-.02604 

-3.1633 

0 

0 

0 

0 

0 

0 

-.142 

-17.25028 

-5.4532 

0 

0 

0 

0 

-.025 

0 

-25.0000 

0 

-.025 

-25.0 

-32.2 

0 

-6.3408 x 10" 4 

0 

-1.205 

5.6506 x 10 2 

- .63408 x 10 4 

-1.205 

0 

1.03178 

0 

0 

0 

0 

0 

0 

0 

1.0 

-2.292 x 10' 3 

0 

-7.0672 

-1.1112 

-2.292 x 10' 3 

-7.0672 

0 

-2.06314 
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TABLE XI 
a 12 MATRIX 

0 

0 

0 

o 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-8.944 x 10' 3 

4.591 x 10" 3 

-2.1262 x 10" 4 

1.5266 x lO -4 

0 

0 

0 

0 

-.1844 

.20762 

-70449 x 10' 3 

6.9711 x 10" 3 



TABLE XI T 


n n 

r\ a 

A 21/u m ' mix 



0 0 o o o O 

o o 0 0 0 0 0 


-1.4343xl0~ 3 

0 

-737.04 

-137.51 

-1.4353xI0" 3 

-738.04 

0 

-133.038 

-3.9012xl0 -3 

0 

752.39 

44.3375 

-3.9017xl0 3 

757.39 

0 

56.4903 
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TABLE XIII 


A 22 ^ matrix 


0 

0 1.0 

0 

0 

0 0 

1.0 

-177.444 

61.8964 -1.159 

.91536 

6.98597 

-456.528 -.121926 

-.849 


TABLE XIV 



B^ MATRIX 


0 0 0 

0 0 -.28852 0 

-15.465 


TABLE XV 



B,', MATRIX 

L, y 


IT - 

CT -2229.4 615.545 



TABLE XVI 



C : MATRIX 



0 0 


0 0 
0 0 0 0 


57.2058 
0 


0 

57.2958 
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TABLE XVII 
CyiATRIX 

-1.4324 =TT6'6158 0 “0 

0 0 -1.4324 -1.66158 


TABLE XVIII 
EJ MATRIX 


5.088 0 0 

0 1 .01948 


0 0 0 0 0 
.10621 OQQ 0 



APPENDIX B 


DERIVATIONS RELATED TO SINGULAR PERTURBATION 


Consider a singularly perturbed linear time- invariant system 


X 1 = A 11 X 1 + A 12 X 2 

+ B x u + EjW, XjCo) = x 1Q 

(A. la) 

wX 2 = A 21 X 1 + A 22 X 2 

+ B 2 u + E 2 w » x 2 (0) = x 20 

(A. lb) 

y = CjXj + C 2 x 2 + 

V 

(A.lc) 


where x^, X 2 , and y are n^, ^ and m dimensional vectors respec- 
tively, the control u is an r vector, and u>0 is a small scalar 
parameter which arises due to the presence of high frequency elastic 
modes. The covariances of the p and q dimensional white noise 
vectors w and v are 


E (w(t) w'(t)} = WS(t-T) 

E {w(t) v “(t) } = 0 
E {v(t) V'(t)) = V<5(t-x) 

Given the observation y(x) for 0<x<~, it is desired to estimate x^(t) 
which is a slowly varying dynamics vector of the system (A.l). 

The conditions that guarantee that high frequency oscillations 
will occur are n^ is even and A ^ has the form [22] 


A 22 


uD. 


1 


D, 


D. 


mD, 


3 

where D 2 , D^ are n ^ /2 x nonsingular matrices and the matrix D 2 D 2 


has simple and negative eigenvalues -u>f, i=l,2,...., n ^ /2. 
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By using the techniques presented in [22-24] it will be shown 
that in the limit (p- + 0 + ) x 2 can be approximated as a white noise 
process which can be used as an input to the slow mode x^. Then the 
estimation of x^(t) by ignoring the high frequency oscillations can 
be analytically represented by the filtering of a reduced order system. 
If n in (A. lb) is neglected and (A. lb) is replaced by 
0 = A 21 x x + A 22 x 2 + B 2 G + E 2 w 

then, if A 22 exists, 

*2 = " A 22 ( A 21 x 1 + B 2 ^ + E 2 W ^ 

and the substitution of x 2 in (A. la) results in the reduced order system 

x, = Ajc, + B G + E w (A. 3a) 

1 o 1 o o v 


y = C x, + D u + F w + v 
o 1 o o 


(A. 3b) 


where 


,-l 


A o A 11 ‘ A 12 A 22 A 21 
•1 


B 0 = B 


E o = E, - 


A 12 A 22 B 2 
A 12 A 22 E 2 


,-l 


C- = C, - A 


21 


D o ~ C 2 A 22 B 2 


F = 


" C 2 A 22 E 2 


without any input, the slowly varying part of x 2 is x 2 = '^22 A 21 X 1' 
separate x 2 from the highly oscillatory part of x 2 , a change of variable 
is used. 

n = x 2 + A 22 A^ x^ + yGx^=x 2 + LXj (A. 4) 
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transforming (A.l) into 

~ ^ A^ 2 n ^ ^ HjW 

un = Fx^ + (A 22 + pLA 12 )n + (B 2 + pLB^p + (E 2 + uLE 1 )w 

where 

F = p(A 22 A 21 + yG) (A q - pA 12 G) - pA 22 G 


Hie solution of F = 0 is 

G = A 22 A 21 A o + °M 
To separate the slow modes, introduce 
C= x^ - y(A^2 A 22 + w ^ n " x i ’ 

and choose M such that 

A 12 + > J (A o - pA 12 G)H - H(A 22 + pLA 12 ) = 0 
so 

M = A q A 12 A 22 - A i2 A 22 A 21 A 12 A 22 + 

The transformation (A. 4) and (A. 5) can be written as 




1 

t 

3 

i 

xr 

F 

-pH" 


X 1 

n 


L 

I 


x~ 



_ 

mj 


L 2J 


The original system (A.l) is finally transformed into 

S = V + V- + V 

pn = A 2 n + B^u + E^w 
y = C Q r, + C 2 n + v 


(A. 5) 


(A. 6) 


(A. 7a) 
(A. 7b) 


(A. 7c) 
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where 

A o = A o ‘ yA 12 G “ A o + 0(p) 

B Q = B q - y (HLB 1 = MB 2 ) « B q - 0(y) 

= E q - U (HLB 1 + ME 2 ) * E q - 0(y) 

A Z = A 22 + yLA 12 ~ A 22 + 

S 2 = B 2 + yLB^ = B 2 + 0(y) 

E 2 - E 2 + yLE^ « E 2 + 000 

= C 1 - C 2 L “ C o + 

C 2 = C 2 + uCq - C 2 )H a C 2 + 0(u) 

To investigate the behavior of n(t), we assume that (A. 7b) can 
be written as 


yn = a 2 n + ff '2 w 


(A. 8) 


That is, u has been replaced by the feedback control before applying 
the transformation to get (A. 7) 
t - t 

Let t = — , then (A* 8) becomes 


d 

ar 


n = <^2 n+ ff 2 w 


where w in the t- scale has covariance W/y instead of W in the t-scale. 
Consequently the covariance of the process n(x) would also have the 
form V/y where V satisfies the eqn. 

H?V = 4 2 V + V^* ¥*2 

with steady state value V satisfying 

0 = i4„V + V A' + EJtE' 

2 » ®> 2 2 2 



If >?2 sta ble and tfivcn an arbitrary c>(), then there exists a*>0 
and tj>0 such that ||V(t) - VJ|<e for all t > tj and ()<a<y*. There- 
fore, for t>tp n(t) may be approximated by a stationary stochastic 
process with the autocorrelation function 

R n (t;r)= e {(O CO) = ^ exp [/i 2 r) ] , r,<r 

However, R (t', tT)->-0 as u->0 for t' ^ t~ , and 

lim/*' 5 R (t' - t~) dr = -A' 1 V - Y/iV 1 
y -»0 t_c n z z 

where 

V w = V + 0(a) 
so that V satisfies 


A 22 V * ^ A 22 + E 2'*2 ■ 0 


It follows that in the limit n becomes a white noise process with 
covariance 


Jio R n (t ' J r) = _(A 22 V + V a 22 5 6(t '‘ r) 


( A 22 E 2 W; 2 A 22 ^ 6 ( t "* O 


which is the covariance of the process 


n A 22 E 2 W ~ A 2 E 2 w 


obtained by the formal substitution of y=0 in (A. 8). 
The estimation equation for the slow mode is 
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Since the inverse of (A. 6) is 



So the filtered estimate of and X2 are given as 
x 1 = c - pHn = t 
x 2 = -L c + (I - pLH)n 

Since /l = A + 0(y), B = B + 0(p),..., so 
o 0 o o 

X 1 = X x + 0(u) 

A 

where x are the estimates obtained by solving the filtering problem for 
the reduced system (A. 3) . 
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